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Abstract 

The c-axis charge dynamics of copper oxide materials in the underdoped and 
optimally doped regimes has been studied by considering the incoherent in- 
terlayer hopping. It is shown that the c-axis charge dynamics for the chain 
copper oxide materials is mainly governed by the scattering from the in-plane 
fluctuation, and the c-axis charge dynamics for the no-chain copper oxide 
materials is dominated by the scattering from the in-plane fluctuation incor- 
porating with the interlayer disorder, which would be suppressed when the 
holon pseudogap opens at low temperatures and lower doping levels, leading 
to the crossovers to the semiconducting-like range in the c-axis resistivity and 

the temperature linear to the nonlinear range in the in-plane resistivity. 
PACS: 71.27.+a, 72.10.-d, 74.72.-h 
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I. INTRODUCTION 



It has become clear in the past several years that copper oxide materials are among 
the most complex systems studied in condensed matter physics, and show many unusual 
normal-state properties. The complications arise mainly from (1) strong anisotropy in the 
properties parallel and perpendicular to the Cu02 planes which are the key structural ele- 
ment in the whole copper oxide superconducting materials, and (2) extreme sensitivity of the 
properties to the compositions (stoichiometry) which control the carrier density in the Cu02 
plane [|T|, while the unusual normal-state feature is then closely related to the fact that these 
copper oxide materials are doped Mott insulators, obtained by chemically adding charge car- 
riers to a strongly correlated antiferromagnetic (AF) insulating state, therefore the physical 
properties of these systems mainly depend on the extent of dopings, and the regimes have 
been classified into the underdoped, optimally doped, and overdoped, respectively 0. The 
normal-state properties of copper oxide materials in the underdoped and optimally doped 
regimes exhibit a number of anomalous properties in the sense that they do not fit in the 
conventional Fermi-liquid theory , and the mechanism for the superconductivity in cop- 
per oxide materials has been widely recognized to be closely associated with the anisotropic 
normal-state properties @J§]. Among the striking features of the normal-state properties in 
the underdoped and optimally doped regimes, the physical quantity which most evidently 
displays the anisotropic property in copper oxide materials is the charge dynamics P] , which 
is manifested by the optical conductivity and resistivity. It has been show from the experi- 
ments that the in-plane charge dynamics is rather universal within the whole copper oxide 
materials The in-plane optical conductivity for the same doping is nearly materials 

independent both in the magnitude and energy dependence, and shows the non-Drude be- 
havior at low energies and anomalous midinfrared band in the charge-transfer gap, while 
the in-plane resistivity Pab{T) exhibits a linear behavior in the temperature in the optimally 
doped regime and a nearly temperature linear dependence with deviations at low temper- 
atures in the underdoped regime By contrast, the magnitude of the c-axis charge 
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dynamics in the underdoped and optimally doped regimes is strongly materials dependent, 
i.e., it is dependent on the species of the building blocks in between the Cu02 planes p|-p!0|. 
In the underdoped and optimally doped regimes, the experimental results p|-p!0| show that 
the ratio R = pdT) / pabiT) ranges from R ~ 100 to i? > 10^, this large magnitude of the 
resistivity anisotropy reflects that the c-axis mean free path is shorter than the interlayer dis- 
tance, and the carriers are tightly confined to the Cu02 planes, and also is the evidence of the 
incoherent charge dynamics in the c-axis direction. For the copper oxide materials without 
the Cu-0 chains in between the Cu02 planes P], such as La2-aSr2.Cu04 systems, the trans- 
ferred weight in the c-axis conductivity forms a band peaked at high energy uj ~ 2ef , and 
the low-energy spectral weight is quite small and spread over a wide energy range instead of 
forming a peak at low energies, in this case the behavior of the c-axis temperature dependent 
resistivity Pc{T) is characterized by a crossover from the high temperature metallic-like to 
the low temperature semiconducting-like [^]. However, for these copper oxide materials with 
the Cu-0 chains in between the Cu02 planes [0, such as YBa2Cu307 systems, the c-axis 
conductivity exhibits the non-Drude behavior at low energies and weak midinfrared band, 
moreover, this weak midinfrared band rapidly decrease with reducing dopings or increasing 
temperatures, while the c-axis resistivity Pc{T) is linear in temperatures in the optimally 
doped regime, and shows a crossover from the high temperature metallic-like behavior to 



the low temperature semiconducting-like behavior in the underdoped regime JTOl . Therefore 
there are some subtle differences between the chain and no-chain copper oxide materials. 
The c-axis charge dynamics of copper oxide materials has been addressed from several 



theoretical viewpoints [jTT]-|l^. Based on the concept of dynamical dephasing, Leggett |TT 
thus proposed that the c-axis conduction has to do with scatterings from in-plane thermal 
fluctuations, and depends on the ratio of the interlayer hopping rate of Cu02 sheets to 
the thermal energy. While the theory of tunneling c-axis conductivity in the incoherent 
regime has been given by many researchers . Based on a highly anisotropic Fermi-liquid, 
some effect from the interlayer static disorder or dynamical one has been discussed |]T^. The 
similar incoherent conductivity in the coupled fermion chains has been in more detail studied 



by many authors within the framework of the non-Fermi-hquid theory |T^. Moreover, the 
most rehable result for the c-axis charge dynamics from the model relevant to copper oxide 



materials has been obtained by the numerical simulation |T5| . It has been argued that the in- 
plane resistivity deviates from the temperature linear behavior and temperature coefficient of 
the c-axis resistivity change sign, showing semiconducting-like behavior at low temperatures 
are associated with the effect of the pseudogap To shed light on this issue, we, in 

this paper, apply the fermion-spin approach [ll6| , p!7| to study the c-axis charge dynamics by 
considering the interlayer coupling. 

The paper is organized as follows. The theoretical framework is presented in Sec. II. In 
the case of the incoherent interlayer hopping, the c-axis current-current correlation function 
(then the c-axis optical conductivity) is calculated in terms of the in-plane single-particle 
spectral function by using standard formalisms for the tunneling in metal-insulator-metal 
junctions Within this theoretical framework, we discuss the c-axis charge dynamics 

of the chain copper oxide materials in Sec. III. It is shown that the c-axis charge dynamics of 
the chain copper oxide materials is mainly governed by the scattering from in-plane charged 
holons due to spinon fluctuations, and the behavior of the c-axis resistivity is the metallic-like 
in the optimally doped regime and the semiconducting-like in the underdoped regime at low 
temperatures. In Sec. IV, the c-axis charge dynamics of the no-chain copper oxide materials 
is discussed. Our result shows that the scattering from the in-plane fluctuation incorporating 
with the interlayer disorder dominates the c-axis charge dynamics for the no-chain copper 
oxide materials. In this case, the c-axis resistivity exhibits the semiconducting- like behavior 
in the underdoped and optimally doped regimes at low temperatures. Sec. V is devoted to 
a summary and discussions. Our results also show that the crossover to the semiconducting- 
like range in Pc{T) is obviously linked with the crossover from the temperature linear to the 
nonlinear range in pab{T), and the common origin for these crossovers is due to the existence 
of the holon pseudogap at low temperatures and lower doping levels. 
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II. THEORETICAL FRAMEWORK 



Among the microscopic models the most simplest for the discussion of doped Mott in- 



sulators is the t-J model |20], which is originally introduced as an effective Hamiltonian of 
the Hubbard model in the strong coupling regime, where the electron become strongly cor- 
related to avoid the double occupancy. The interest in the t-J model is stimulated by many 
researchers' suggestions that it may contain the essential physics of copper oxide materials 
2^,^. On the other hand, there is a lot of evidence from the experiments and numerical 



simulations in favour of the t-J model as the basic underlying microscopic model P,p2|. 
Within each Cu02 plane, the physics property may be described by the two-dimensional 
(2D) t-J model. 

Hi = -t^ Cl^Cli+fja- + h.C. - /i X! ^lia^lio- + J^^H' S^j+j), (l) 
ifja icr irj 

supplemented by the on-site local constraint J^a Cli^Cuu < 1 to avoid the double occupancy, 
where f] = ±aoX, ^aoV, clq is the lattice constant of the square planar lattice, which is set 
as the unit hereafter, i refers to planar sites within the 1-th Cu02 plane, C^^ {Cua) are 
the electron creation (annihilation) operators, Su = C/jCrC;i/2 are the spin operators with 
cr = {ax,<Jy,<7z) as the Pauh matrices, and fj. is the chemical potential. Then the hopping 
between Cu02 planes is considered as [TH 



H = -tcY: CLCi+f,.^. + h.C. + Y: Hi, (2) 

Irjcicr I 

where = icgi, Cq is the interlayer distance, and has been determined from the experiments 
p3| as Co > 2ao. As mentioned above, the experimental results show that the c-axis charge 



dynamics in the underdoped and optimally doped regimes is incoherent, therefore the c-axis 
momentum can not be defined ||2^. Moreover, the absence of the coherent c-axis charge 
dynamics is a consequence of the weak interlayer hopping matrix element tc, but also of a 
strong intralayer scattering, i.e., tc -C t, and therefore the common Cu02 planes in copper 
oxide materials clearly dominate the most normal-state properties. In this case, the most 



5 



relevant for the study of the c-axis charge dynamics is the results on the in-plane conductivity 
cTafe(co') and related single-particle spectral function A{k,u). 

Since the strong electron correlation in the t-J model manifests itself by the electron sin- 
gle occupancy on-site local constraint, then the crucial requirement is to impose this electron 
on-site local constraint for a proper understanding of the physics of copper oxide materials. 
To incorporate this local constraint, the fermion-spin theory based on the charge-spin sep- 
aration has been proposed |[16| , P^ . According to the fermion-spin theory, the constrained 



electron operators in the t-J model is decomposed as ||16 



Clii - ^li^ii ' ^lil - ^li^li^ (3) 

with the spinless fermion operator hi keeps track of the charge (holon) , while the pseudospin 
operator Si keeps track of the spin (spinon). The main advantage of this approach is that 
the electron on-site local constraint can be treated exactly in analytical calculations. In this 
case, the low-energy behavior of the t-J model (2) in the fermion-spin representation can be 
written as [ ^ , 

H = t,Y. hU^MiSt.S^^^^, + Sr^St^^J + E Hu (4a) 

Ifj^i I 

Hi = tY^ hii^fjhii{Sj^ Sfi_^fj + Sj^S^^fj) + fJ'Yl ^h^ii + "^e// ■ ^^+'7)' (4b) 

ifi i if) 

where Jeff = J[{1 — — cj)'^], the in-plane holon particle-hole parameter cj) = {hl^hu^fi), and 
and Sj^ are the pseudospin raising and lowering operators, respectively. As a consequence, 
the kinetic part in the t-J model has been expressed as the holon-spinon interaction in 
the fermion-spin representation, which dominates the charge and spin dynamics in copper 
oxide materials in the underdoped and optimally doped regimes. The spinon and holon 
may be separated at the mean-field level, but they are strongly coupled beyond mean-field 
approximation (MFA) due to fluctuations. 

The mean-field theory within the fermion-spin formalism in the underdoped and op- 
timally doped regimes without AF long-range-order (AFLRO) has been developed []T7 
and the in-plane mean-field spinon and holon Green's functions D^abi^ ~ J? ^ — t') = 



— {TrS^{T)Sij{T'))o and g'^h{i — j, t — t') = —{Trhii{T)hjj{T'))(j have been evaluated |T7 
as, 

respectively, where Bk = X[{2eXz + x)lk - (ex + 2%^)], 7k = (1/^) E-^ e*^'', = '^ZJ^ff, 
e = 1 + 2t(f)/Jf;ff, Z is the number of the nearest neighbor sites at the plane, while the 
in-plane mean-field spinon spectrum 

u\k) = \^ (ae[xzjk + ^X] - [(^C, + ^(1 - «)]) i^lk - 1) 

+X' (aei^xik + ^Xz] - l^[aC + ^(1 - a)]) (7^ - e), (7) 

and the in-plane mean-field holon spectrum = '^Zxt'fk + A^, with the in-plane spinon 
correlation functions x = {ShSh+^), Xz = {SuSti+i^), C = (1/^^) E,),^"/ ('^^j+^S^T^^',), and 
Cz = iX/Z"^) J2fj ^1 {Sii+fi^ii_^^i) ■ order not to violate the sum rule of the correlation function 
{Sj^Sn) = 1/2 in the case without AFLRO, the important decoupling parameter a has 



been introduced in the mean- field calculation [I^, which can be regarded as the vertex 
correction. The mean-field order parameter x, C, Xz, Cz, and chemical potential n have 
been determined ||T3 by the self-consistent equations. 

Within the 2D t-J model, the in-plane charge dynamics in copper oxide materials has 
been discussed ||2^ by considering fiuctuations around this mean-field solution, and the result 
exhibits a behavior similar to that seen in the experiments 0,^] and numerical simulations 
p6| . In the framework of the charge-spin separation, an electron is represented by the 
product of a holon and a spinon, then the external field can only be coupled to one of them. 
According to the loffe-Larkin combination rule P7| , PS| , the physical c-axis conductivity (Jc{uj) 
is given by, 

a;\u;) = ai^^'\u) + ai^^~\u), (8) 



where a^'^^uj) and a^'^\uj) are the contributions to the c-axis conductivity from holons and 
spinous, respectively, and can be expressed JTSi as. 



with and {uj) are the holon and spinon c-axis current- current correlation function, 

respectively, which are defined as, 

nW(r - r') = -{Trji'\r)j^'\T')), nW(r - r') = -(T.jf Mjf )). (10) 

Within the Hamiltonian (4), the c-axis current densities of spinous and holons are obtained 
by the time derivation of the polarization operator using Heisenberg's equation of motion 

as, jf) = 2icexEif,,iVchl+fj^ihii and j^'^ = tce(f)cEi^,iVc{Su Sf^^^^j^i + 5^7^^+^^^, respectively, 
with tc = tcXc/x is the effective interlayer holon hopping matrix element, and the order 
parameters Xc and 0c are defined fl^ as Xc = (SuSi^fi^), and 0c = (hlihi^fjci) y respectively. 



As in the previous discussions p5| , a formal calculation for the spinon part shows that there 
is no the direct contribution to the current-current correlation from spinous, but the strongly 
correlation between holons and spinous is considered through the spinon's order parameters 
entering in the holon part of the contribution to the current-current correlation, therefore 
the charge dynamics in copper oxide materials is mainly caused by charged holons within 
the Cu02 planes, which are strongly renormalized because of the strong interaction with 
fluctuations of surrounding spinon excitations. 

In the case of the incoherent charge dynamics in the c-axis direction, i.e., the independent 
electron propagation in each layer, the c-axis holon current-current correlation function is 
then proportional to the tunneling rate between just two adjacent planes, and can be calcu- 
lated in terms of the in-plane holon Green's function gab{k, uj) by using standard formalisms 
for the tunneling in metal-insulator- metal junctions |[T5| , p!8|J19[| as, 

Uf\iujn) = -{dcexcof^Y.^Y.9ab{k,iuj'^ + iuJn)gab{k,iu'^), (11) 

k P ico',^ 

where iun is the Matsubara frequency. Therefore the c-axis current-current correlation 
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function is essentially determined by the property of the in-plane full holon propagator 
gab{k,iUn), which can be expressed as, 

gab~^{k,iuJn) -T.^h\k,iuJn) iuJn - - ^h\k, iuJnY ^ ^ 

where the holon self-energy has been obtained by considering the second-order correction 
due to the spinon pair bubble as [p5| , 

( c^nF{ip+k)[nB{^p') - nB{uJp+p')] - nB{u}p+p')nB{-ujp') 



iUn + UJp+p' — UJp' — ^p+k 

nF{^p+k)[nB{i^p+p') -nB{-uJp')] + nB{uJp')nB{uJp+p' 



iuJn + iOp/ + UJp-i-p' — ^p+k 

nF{^p+k)[nB{i^p+p') - nB{-uJp:)] - nB{-Up')nB{-u]p+p' 



(13) 



ioJn — UJp+p' — i^p' — ^p+k J 

with nF{^k) and nB{uJk) are the Fermi and Bose distribution functions, respectively. For 
the convenience in the following discussions, the above full holon in-plane Green's function 
Qabik, iojn) also can be expressed as frequency integrals in terms of the spectral representation 



gab{k,iujn)= 7^ — : , (14) 

J-oo /Vr tUJn — 00 

with the in-plane holon spectral function Afi{k,uj) = --2lm.gab{k,uj). Then the c-axis optical 
conductivity in the present theoretical framework is expressed [13] as aduj) = 



III. CHARGE DYNAMICS FOR THE CHAIN COPPER OXIDE MATERIALS 

We firstly consider the chain copper oxide materials. From the experiments testing the 
c-axis charge dynamics [l^, it has been shown that the presence of the rather conductive 
Cu-0 chains in the underdoped and optimally doped regimes can reduce the blocking effect, 
and therefore the c-axis charge dynamics in this system is effected by the same electron 
interaction as that in the in-plane. In this case, we substitute Eq. (14) into Eq. (11), and 
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evaluate the frequency summations, then the c-axis optical conductivity for the chain copper 
oxide materials can be obtained as p!5|JT8|JT9[] , 

= -(4t,exco)'-E / ^Ah{k,J^uj)Ah{k,J)^ (15) 

where the in-plane momentum is conserved. This c-axis conductivity (Jc{uj) has been calcu- 
lated numerically, and the result at the doping 6 = 0.12 (solid line), 6 = 0.09 (dashed line), 
and 6 = 0.06 (dot-dashed line) for the parameters t/J = 2.5, tc/t = 0.04, and co/cto = 2.5 
at the temperature T=0 is shown in Fig. 1, where the charge e has been set as the unit. 
From Fig. 1, we find that there is two bands in aduj) separated at ~ 0.4t, the higher- 
energy band, corresponding to the "midinfrared band" in the in-plane optical conductivity 
aab{^^) P,|3|,|25[| , shows a broad peak at c<j ~ 0.7t, in particular, the weight of this band is 
strongly doping dependent, and decreasing rapidly with dopings, but the peak position does 
not appreciably shift to higher energies. On the other hand, the transferred weight of the 
lower-energy band forms a sharp peak at u; < OAt, which can be described formally by 
the non-Drude formula, and our analysis indicates that this peak decay is l/u at low 
energies as in the case of <7ab{^^) P,|3|J25[]. In comparison with o"afe(c<j) |2^, the present result 
also shows that the values of (Jc{uj) are by 2 ~ 3 orders of magnitude smaller than those of 
aab{^^) in the corresponding energy range. The finite temperature behavior of aduj) also has 
been discussed, and the result shows that crc(^) is temperature dependent, the higher-energy 
band is severely suppressed with increasing temperatures, and vanishes at higher temper- 
atures. These results are qualitatively consistent with the experimental results [0] of the 
chain copper oxide materials and numerical simulations |1^ . 



With the help of the c-axis conductivity, the c-axis resistivity can be obtained as 
Pc = 1/ lim^_+o crc('^)5 and the numerical result at the doping 6 = 0.12 and 6 = 0.06 for 
the parameters t/J = 2.5, ic/t = 0.04, and Co/ao = 2.5 is shown in Fig. 2(a) and Fig. 
2(b), respectively. In the underdoped regime, the behavior of the temperature dependence 
of Pc{T) shows a crossover from the high temperature metallic-like {dpc{T)/dT > 0) to the 
low temperature semiconducting- like {dpc(T)/dT < 0), but the metallic-like temperature 
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dependence dominates over a wide temperature range. Therefore in this case, there is a gen- 
eral trend that the chain copper oxide materials show nonmetallic Pc{T) in the underdoped 
regime at low temperatures. While in the optimally doped regime, pc{T) is linear in tem- 
peratures, and shows the metallic-like behavior for all temperatures. These results are also 
qualitatively consistent with the experimental results of the chain copper oxide materials 
||T0| and numerical simulation |^ . 



IV. CHARGE DYNAMICS FOR THE NO-CHAIN COPPER OXIDE MATERIALS 

Now we turn to discuss the c-axis charge dynamics of the no-chain copper oxide materials. 
It has been indicated from the experiments that for the no-chain copper oxide materials 
the doped holes may introduce a disorder in between the Cu02 planes, contrary to the case 



of the chain copper oxide materials |T0[ , where the increasing doping reduces the disorder in 
between the Cu02 planes due to the effect of the Cu-0 chains. Therefore for the no-chain 
copper oxide materials, the disorder introduced by doped holes residing between the Cu02 
planes modifies the interlayer hopping elements as the random matrix elements. In this 
case, only the in-plane holon density of states (DOS) flh{^) = '^/NJ2k^h{k,u!) enters the 
holon current-current correlation function as in disordered systems ||T5| , ^ , and after a similar 
discussion as in Sec. II, we find that the corresponding momentum-nonconserving expression 
of the c-axis conductivity cr^"^ (to) for the no-chain copper oxide materials is obtained by the 
replacement of the in-plane holon spectral function Ah{k,uj) in Eq. (15) with the in-plane 
holon DOS nhicu) as Pj2|], 



<yT\^) = o(4tcexco)' / —VthiJ + uj)Vth{uj')^ '- (16) 

where the if. is some average of the random interlayer hopping matrix elements {t^u. We 
have performed a numerical calculation for this c-axis conductivity cr^")(co'), and the result 
at the doping 5 = 0.10 (solid line), 6 = 0.08 (dashed line), and 6 = 0.06 (dash-dotted line) 
for the parameters t/J = 2.5, ic/t = 0.04 with the temperature T = is plotted in Fig. 3. 
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This result shows that the c-axis conductivity 0"^"^ (cj) contains two bands, the higher-energy 
band, corresponding to the midinfrared band in the in-plane conductivity aab{i^) [Pi|3|,|25| , 
shows a broad peak at ~ 0.3t. The weight of this band is increased with dopings, but the 
peak position does not appreciably shift to lower energies. As a consequence of this pinning 
of the transferred spectral weight, the weight of the lower-energy band, corresponding to 
the non-Drude peak in aab{uj) P,p|,|25| , is quite small and does not form a well-defined 
peak at low energies in the underdoped and optimally doped regimes. In this case, the 
conductivity a^^^to) at low energies cannot be described by the over-damped Drude-hke 
formula even an w-dependence of the in-plane holon scattering rate has been taken into 
consideration. In comparison with the momentum-conserving crd^), our result also shows 
that at low energies the suppression of the momentum-nonconserving a^"-\uj) is due to 
the disordered effect introduced by doped holes residing between the Cu02 planes. These 
results are qualitatively consistent with the experimental results of the no-chain copper oxide 
materials 0. 

For the further understanding the transport property of the no-chain copper oxide 
materials, we have also performed the numerical calculation for the c-axis resistivity 
p^"-* = 1/ lim^_^o Cc^H"^)) the results at the doping S = 0.10 and 5 = 0.06 for the 
parameters t/J = 2.5, tc/t = 0.04, and Cq = 2.5ao are shown in Fig. 4(a) and Fig. 4(b), 
respectively. In accordance with the c-axis conductivity a^"'\uj), the behavior of the c-axis 
resistivity p^^\T) in the underdoped and optimally doped regimes is the semiconducting- 
like at low temperatures, and metallic-like at higher temperatures. In comparison with the 
in-plane resistivity Pab{T) it is shown that the values of the c-axis resistivity p^^\T) for 
the no-chain copper oxide materials are by 3 ~ 4 orders of magnitude larger than these of 
the in-plane resistivity Pab{T) in the corresponding energy range, which are also qualitatively 
consistent with the experimental results of the no-chain copper oxide materials [^. 
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V. SUMMARY AND DISCUSSIONS 



In the above discussions, the central concerns of the c-axis charge dynamics in copper 
oxide materials are the two dimensionality of the electron state and incoherent hopping be- 
tween the Cu02 planes, and therefore the c-axis charge dynamics in the present fermion-spin 
picture based on the charge-spin separation is mainly determined by the in-plane charged 
holon fluctuation for the chain copper oxide materials and the in-plane charged holon fluc- 
tuation incorporating with the interlayer disorder for the no-chain copper oxide materials. 
In comparison with the in-plane resistivity Pab{T) it is shown that the crossover to the 
semiconducting-like range in Pc{T) is obviously linked with the crossover from the temper- 
ature linear to the nonlinear range in pab{T), i.e., they should have a common origin. 

In the fermion-spin theory [0, the charge and spin degrees of freedom of the physical 
electron are separated as the holon and spinon, respectively. Although both holons and 
spinous contributed to the charge and spin dynamics, it has been shown that the scattering 
of spinous dominates the spin dynamics [|30|, while the results of the in-plane charge dynamics 
p5| and present c-axis charge dynamics show that scattering of holons dominates the charge 
dynamics, the two rates observed in the experiments [0] are attributed to the scattering of 
two distinct excitations, spinous and holons. It has been shown that an remarkable point 
of the pseudogap is that it appears in both of spinon and holon excitations [P]-pi]|]. The 
present study indicates that the observed crossovers of pab and pc for copper oxide materials 
seem to be connected with the pseudogap in the in-plane charge holon excitations, which 
can be understood from the physical property of the in-plane holon DOS. The numerical 
result of the in-plane holon DOS ^lh{^) at the doping 6 = 0.06, 6 = 0.12, and 6 = 0.15 
for the parameter t/J = 2.5 at the temperature T=0 is shown in Fig. 5(a), Fig. 5(b), and 
Fig. 5(c), respectively. For comparison, the corresponding mean- field result (dashed line) is 
also shown in Fig. 5. While the in-plane holon density of states ^hi^^) in the underdoping 
6 = 0.06 as a function of energy for the temperature (a) T = 0, (b) T = O.OIJ, and (c) 
T = 0.2J is plotted in Fig. 6. From Fig. 5 and Fig. 6, we therefore find that the in-plane 
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holon DOS in MFA consists of the central part only, which comes from the noninteracting 
particles as pointed in Ref. |2^. After including fluctuations the central part is renormalized 
and two side bands |^ and a V-shape holon pseudogap near the chemical potential /i in the 



underdoped regime appear. But these two side bands are almost doping and temperature 
independent, while the V-shape holon pseudogap is doping and temperature dependent, 
and grows monotonously as the doping 6 decreases, and disappear in the overdoped regime. 
Moreover, this holon pseudogap also decreases with increasing temperatures, and vanishes 
at higher temperatures. Since the full holon Green's function (then the holon spectral 
function and DOS) is obtained by considering the second-order correction due to the spinon 
pair bubble, then the holon pseudogap is closely related to the spinon fluctuation. For small 
dopings and lower temperatures, the holon kinetic energy is much smaller than the magnetic 
energy, i.e., 6t ^ J, in this case the magnetic fluctuation is strong enough to lead to the 
holon pseudogap. This holon pseudogap would reduce the in- plane holon scattering and thus 
is responsible for the metallic to semiconducting crossover in the c-axis resistivity pc and the 
deviation from the temperature linear behavior in the in-plane resistivity pab This holon 



pseudogap will also lead to form the normal-state gap in the system, and the similar result 



has been obtained from the doped kagome and triangular antiferromagnets PTI], where the 
strong quantum fluctuation of spinous due to the geometric frustration leads to the normal- 
state gap. With increasing temperatures or dopings, the holon kinetic energy is increased, 
while the spinon magnetic energy is decreased. In the region where the holon pseudogap 
closes, at high temperatures or at higher doping levels, the charged holon scattering would 
give rise to the temperature linear in-plane resistivity as well as the metallic temperature 
dependence of the c-axis resistivity. Our results also show that PabiT) is only slightly affected 
by this holon pseudogap [^, while pc{T) is more sensitive to the underlying mechanism. 

In summary, we have studied the c-axis charge dynamics of copper oxide materials in 
the underdoped and optimally doped regimes within the t-J model by considering the in- 
coherent interlayer hopping. Our result shows the c-axis charge dynamics for the chain 
copper oxide materials is mainly governed by the scattering from the in-plane fluctuation, 
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and the c-axis charge dynamics for the no-chain copper oxide materials is dominated by the 
scattering from the in-plane fluctuation incorporating with the interlayer disorder, which 
would be suppressed when the holon pseudogap opens at low temperatures and lower dop- 
ing levels, leading to the crossovers to the semiconducting-like range in the c-axis resistivity 
Pc(T) and the temperature linear to the nonlinear range in the in-plane resistivity pab{T). 
Because copper oxide materials are very complex systems, it is also possible that the actual 
c-axis conductivity may be a linear combination of the momentum-conserving adcu) and 
momentum-nonconserving a^^\uj), but we believe that aduj) should be the major part of 
the c-axis conductivity in the chain copper oxide materials, while a^""^ should be the major 
part of the c-axis conductivity in the no-chain copper oxide materials. 
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FIGURES 

FIG. 1. The c-axis optical conductivity for the chain copper oxide materials at the doping 
(5 = 0.12 (soUd line), 8 = 0.09 (dashed line), and 5 = 0.06 (dot-dashed line) for t/J = 2.5, 
ic/t = 0.04, and co/ao = 2.5 with the temperature T = 0. 

FIG. 2. The c-axis electron resistivity pc for the chain copper oxide materials at t/J = 2.5, 
ic/t = 0.04, and co/oq = 2.5 for (a) the doping (5 = 0.12 and (b) S = 0.06. 

(n) 

FIG. 3. The c-axis optical conductivity Uc for the no-chain copper oxide materials at the 
doping 5 = 0.10 (solid line), 6 = 0.08 (dashed line), and 6 = 0.06 (dot-dashed line) for t/J = 2.5, 
ic/t = 0.04, and co/ao = 2.5 with the temperature T = 0. 

FIG. 4. The c-axis electron resistivity pc"^ for the no-chain copper oxide materials at t/J = 2.5, 
ic/t = 0.04, and co/ao = 2.5 for (a) the doping 5 = 0.10 and (b) 5 = 0.06. 

FIG. 5. The in-plane holon density of states at t/J = 2.5 for (a) the doping 6 = 0.06, (b) 
5 = 0.12, and (c) 5 = 0.15. The dashed line is the result at the mean-field level. 

FIG. 6. The in-plane holon density of states for t/J = 2.5 at the doping 6 = 0.06 for the 
temperature (a) T = 0, (b) T = O.OIJ, and (c) T = 0.2J. 
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